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Abstract 

It is shown that the potential for lepton-antilepton bound states (leptonium) 
is the Fourier transform of the first Born approximation to the QED scattering 
amplitude in an 8-component equation, while 16-component equations are 
excluded. The Fourier transform is exact at all cms energies — oo < E < oo; 
the resulting atomic spectrum is explicitly CPT- invariant. 
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I. INTRODUCTION 



In muonium e~/i + , the mass m 2 of the heavier particle is so large that the kinetic energy 
p 2 /2m 2 can be added as a recoil correction to the electron's Dirac Hamiltonian. The resulting 
Dirac equation has 8 components, with the Pauli matrices er 2 in the hyperfine operator . 
However, a similar 8-component equation exists which avoids the expansion in terms of l/m 2 
HUH]. Despite its asymmetric form, it reproduces the energy levels also for positronium e~e + 
to order a A [0]. Unfortunately, its complicated derivation via Breit operators introduced 
a low-energy approximation, which is removed in the present paper. The new interaction 
is simply the Fourier transform of the QED Born term in the cms, valid at all energies, 
— oo < E < oo (by crossing symmetry, E < for e~/i + -scattering refers to E > for 
e + /i~-scattering). The equation is written in the cms, p = Pi = — P2, in units h = c = 1: 

l -E 2 ^ = (K + KjW , K = Epa + l -{m 2 + (3 mi ) 2 , (1) 
« = 750-1 , 75 = ( ? n ) , P = I n _? J ■ (2) 



The Pauli matrices a 2 appear only in the interaction operator Kj. The proof that the 
equation with Kj = describes two free Dirac particles instead of one is given in section [0]. 
The construction of Kj from the Born approximation is done in section |TTT| . Orthogonality 



relations are derived in section IV 



Before going into details, we remind the reader that in a 16-component version, 

Eip = {H x + H 2 + H n )ip , Hi = Pl a, + , (3) 

the energy levels to order a 4 are not reproduced when H\ 2 is the Fourier transform of 
the Born approximation, H\ 2 = V(l — 0*10:2) (V = —a/r, a = e 2 ). This interaction 
was improved by Breit || and was later amended by energy-transfer |^,|| and by positive- 
energy projectors ||. The method is now quite successful in atomic theory, where the 
nucleus produces an external potential for the electrons ||. In the previous derivation of 
(0) via Breit operators a simpler prescription was used which sets the (large) squares 
of Breit operators equal to zero. However, the connection with perturbative QED was then 
incomplete. In the present derivation, the connection with the Born series is straightforward. 
With loop graphs included, all a 6 -effects in binding energies should appear, for arbitrary 
values of m\jm 2 . 

For V = 0, fll]) and (^) are equivalent. In the first Born approximation, Kj of ([]]) and H\ 2 
of AH) are linear in V. When @ is reduced to 8 components, quadratic terms appear which 
contradict the leptonium energy levels to order a 4 . Some of these terms were eliminated by 
a canonical transformation, the rest by projectors. 

To understand the significance of the 8-component theory, one may note that 8 is the 
minimum number of components for a Lorentz- and parity invariant theory. An irreducible 
representation of the Lorentz group for a single lepton requires only two components, for 
example the right-handed i/)r in a chiral basis. It is the parity transformation, r) = 
f3ip(r), which exchanges ipR with the left-handed components tpi (in the chiral or Weyl basis, 
75 and P of (0) exchange their roles). For two leptons, the necessary 2x2 = 4 components 
are doubled by the matrix /3 of (pj), which corresponds to P\P 2 in a 16-component formalism. 
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The separate matrices (3\ and (3 2 of the 16-component formalism provide another unnecessary 
doubling. 

We end this section by casting ([1]) into a more familiar form. We divide it by E and 
define the reduced mass u and reduced energy e: 

u = mi m 2 /E , e={E 2 -m\- mj) /2E , (4) 
ei) = (pa +/1/3 + Kj/E)^ . (5) 



We also anticipate the result (|4l]) for Ki, which gives: 



Kj/E — V — l^Vp((T 1 X <T 2 )/E . (6) 

The difference to an effective Dirac equation is obviously restricted to the hyperfine operator 
— iVp(cr 1 x er 2 )/E. As the cms energy E is close to mi + m 2 for ordinary bound states, the 
replacement of the familiar l/m 2 by l/E is a recoil correction. The essential new feature is 
that the operator Vp is not hermitian. An 8-component equation with a hermitian operator 
does not reproduce the positronium bound states. 

Anomalous magnetic moments are neglected for the time being, so that the present 
formulation does not cover atomic hydrogen. The annihilation graph in the 3 Si state of 
positronium is also neglected. The hyperfine operator will be further simplified in appendix 
[A[ Its connection with the previous approximate operator j|,[5| is given in appendix ||. 
Appendix contains a spin summation, which is included for demonstration of covariance. 



II. EIGHT-COMPONENT SPINORS FOR TWO LEPTONS 

Let Uig and Uif denote the large and small components of the Dirac spinor u of a particle 
i in an orbital: 

PiO-iUig = (-nf + TOi) u if , PitTiUif = (7r° - mi) u ig , 7if = Ei - V e (ri) . (7) 

Here V e is a possible external potential, as in He-atoms. Leptonium has V e = 0, of course. 
The chiral operator 7 5i exchanges u ig with Uif. Out of the 4 products Ui g u 2g ■ ■ ■ Uiju 2 j, we 
keep in this section only the two combinations of total chirality +1 (the eigenvalue of 751752): 

~^{UlgU 2 g + U lf U 2f ) = U d g , 

■^{ui f u 2g + u lg u 2} ) = u df . (8) 

Inspection shows that these combinations satisfy the following equations: 
(m 2 p 1 cr 1 + m x p 2 a 2 ) u dg = (m 2 7r° + m^) u df , 

(m 2 p 1 (T 1 - m x p 2 cr 2 ) u df = (m 2 7r° - m^) u dg . (9) 

Apart from the mi-terms, these equations have the same structure as a single-particle Dirac 
equation. We therefore introduce an 8 x 8-matrix (3 (which has the same eigenvalues as 
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/?i/3 2 ) and an anticommuting matrix 75, 75/? = —fljs, and combine u dg and Udf into one 
8-component double spinor u d : 




u d = ( * ' 75(^2Picn + A™ip 2 cr 2 ) w d = (m 2 7T? - fimiiil) u d . (10) 



75 and /3 have already been given in (|2|), understanding that the number 1 is a 4 x 4- matrix in 
spin space. For mi = 777.2, (HO) should also describe two non-interacting electrons in a helium 
atom. Their Born scattering amplitude for V e 7^ will be constructed in (|35| ) below, but 
our construction of the interaction operator in the differential equation neglects V e . Turning 
now to the free two-body problem, we set pi = — p 2 = p , 

n " = E ' = 2^ 2 + ™?-™2), A = E 2 = ^{E 2 -ml + ml). (11) 
A factor m 2 — (3mi can then be divided off from the right-hand side of (|T0"D, the result being: 

(a c p + fi/3 - e) u d = , (12) 
ot' = j 5 crl = 7 5 (m 2 cri - /3mi<r 2 )/(m 2 + (3mi) . (13) 

Including the space-dependence, the free double-spinor wave function in the cms is: 

^(r) = u d e« kiri e ik2r2 = u d e* kr , r = n - r 2 . (14) 

It differs from the corresponding free-particle solution of ([I]), which has ct c replaced by a. 
This difference may be one reason why (|l|) has not been discovered 60 years ago. During these 
60 years, many different formalisms have been developed. Bethe and Salpeter advocated 
the use of four-dimensional integral equations (with a relative time as fourth integration 
variable), again with 16 components. Although the relative-time concept turned out to be 
useless, one learned to find bound states from integral equations in momentum space, now 
in three dimensions 0. This method avoids the Fourier transformation. Having elaborated 
these momentum space methods, one may be unwilling to return to differential equations, 
particularly if these require such strange matrices as cr\. On the other hand, the advantage 
of the present formulation survives also in 8-component integral equations in momentum 
space. It could become essential also in quarkonium models. 

To establish the connection between (|I|) and (|T2"D, we first define two mass operators: 

m± = m 2 ± f3m 1 , m + m_ = m 2 — ml , m\ = m\ + m 2 ± 2mim 2 /3 . (15) 

The expression for m 2 + has already been used in (|5|). From 75/? = — ^75, one finds: 

7 5 m+ = m_7 5 . (16) 

The Dirac spin operators are 75 cri and 7scr 2 . With the algebra (|15"D, (|T6|), one easily verifies: 

(a c ) 2 = in^crl) 2 = {m 2 cr l + ( 5micr 2 )(m 2 cr 1 - (3m 1 a- 2 )/m + m_ = 3 . (17) 

Comparing this with a 2 = 3, one sees that there should exist a transformation from of to 
a. It was first found in an explicit decomposition of u d into u dg and u d f, and of the spin 
states X12 into \s (singlet) and Xt (triplet, eq. (5.14) in ||). Its compact Dirac form is: 
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ipd = cijj , c = (m + m_) x l 2 [m + — 2mi[3P s ] . (18) 

P s is the projector on the singlet spin state Xs- in the following, the corresponding triplet 
projectors on the three states Xt will also be needed: 

F s = -(l-o^), F t = J(3 + o-i<7 2 ). (19) 
We also define combinations of Pauli matrices, 

o = Oi + o 2 , Act = 01 - a 2 , ct x = ct x x ct 2 , (20) 
which have the following products with P s and F t : 

ctF s = F s ct = , A(TP t = F s Act , ActF s = F t Ao , ct x F s = P t cr x . (21) 
c is not unitary, its inverse being: 

(T 1 = (m+m_)" 1/2 (m + - 2m 1 f3P t ) , (22) 
which is checked by using P s + P t — 1, F s F t = 0. An important property of c is: 

c~ 1 75 = 75C (23) 

By parity invariance, the operators ( p0|) are always accompanied by one factor 75, which in 
view of (p3|) replaces c _1 by c: 

ca c = er m + /m_ , cActc = Act, cct x c = ct x . (24) 

Expressing o - ! as |(ct + Act), one finds for the combination dl3|): 

cct^ c = CTi . (25) 

Thus the transformation (|T8| ) leads from fli~2|) to the free equation (|l]). Other forms of that 
equation are generated by additional transformations d satisfying: 

(3d = d(3 d- 1 75 = 75 rf. (26) 



III. THE BORN APPROXIMATION AND ITS FOURIER TRANSFORM 

The Lorentz-invariant T- matrix for lepton-antilepton scattering from initial orbitals 1,2 
into final orbitals 1', 2' has the following Born approximation: 

T/4tt = au' 1 r y^uiu 2 r y f _ l U2/t = aw'/w'^ (1 — ctia. 2 )uiU2/t , (27) 

with 7^ = /3jO!j, and t = q^q^ = q$ — q 2 being the square of the 4-momentum transfer. In 
the cms, the arguments of the free Dirac spinors are k and — k in the initial state and k' 
and — k' in the final state, and q = k — k' , q = . Remembering 
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mJ uiujtia = {u'lgUig + u'l f uif) (u 2g u 2g + uJ/W2/) (28) 

etc., one sees that T cannot be written as a bilinear in Ud and w^f. One also needs the states 
of total chirality —1, which will be called Wd g and Wdf' 

W dg = -f=(UlgU 2 g - U lf U 2 f) , 

W df = -j={ui}U 2 g - U lg U 2f ) , (29) 
T/AlX = a u'd(l — <Tx<T 2 )Ud + w' d \(1 + <Tx(T 2 )Wd /t. (30) 

This form still has its 16-component character, as the replacement of uiu 2 by Ud and Wd is 
just a unitary transformation. However, in addition to the separate equations for Ud and 
Wd (the equation for Wd has cr 2 replaced by — cr 2 ) there exist also coupled equations, with 
7T° = 7r° + ix 2 and p± = Pio - ! ± p 2 cr 2 : 

w d = m+^vr - -f 5 P+)u d , u d = m^vr - ^p-)w d , (31) 

which can be verified explicitly from (|j) and (|29|). By means of ([31]), T can be written in 
terms of a single 8 x 8-matrix M, 

T/Att = aw' ] d Mud/t, (32) 
M = (tt° - 75 p , _)m; 1 (l - <x lC r 2 ) + (1 + tr^m^ir - l5 p + ) . (33) 

The operators proportional to tt° combine into 27T°, while the operators containing p'_ com- 
bine as follows: 

P '_ (1 - (T \fJ 2 ) = p(A<X + 2<X X ) = p_(l - (Tjffj) . (34) 

The total momentum p = pi + P2 = Pi + P 2 commutes with V and vanishes in the cms. 
There, M reduces to: 

M = m; 1 [2E - 75(1 + cr 1 <T 2 )kAcr] = 1m^{E - z 75 kcr x ) . (35) 
An equivalent form of T follows from the elimination of w'J and u d in (PD|), 

T/4vr = a u'j M f , (36) 
M f = [2£ - 7 5 k'Acr(l + cr lC r 2 )] m^ 1 = 2(E + <r x )m+ 1 . (37) 

This suggests the definition of a second double-spinor wave-function as follows: 

Xd{?) = ^ r m + w d . (38) 



When T is expressed in terms of m + Wd and w'Jm + , the factor m + 1 vanishes both in ( pop 
and (|37|) . The interaction in coordinate space follows as: 

-4tt a/q 2 = / d 3 re-* k ' r V e ikr , -4tt a k/q 2 = J d 3 re~ ik ' r Vp e ikr , (39) 

-4vr a k'/q 2 = J rf 3 re" ik ' r e ikr . (40) 
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It produces the operator 

Kj = V(E - w x 75 ) , K\ = (E + ip<T x l5 )V, (41) 
to be used in (|I|) and in the corresponding equation for the wave-function x(r) = e* kr u>: 

l -E 2 X = (K + K\) X . (42) 

Although Kj is not hermitian, Kj and K\ give equivalent equations, such that the bound 
state energies may be real. 

IV. ORTHOGONALITY RELATIONS AND CONCLUDING REMARKS 

When E is replaced by m = mi + iri2 in the hyperfine operator, (||) and the corresponding 
equation fl42|) assume Hamiltonian forms: 

e^ = H^, e X = H\. (43) 

Taking the hermitian adjoint of the second equation at reduced energy e' and integrating 
over r, one obtains 

(e - e') J X 'H = J xHH -H)i/> = 0. (44) 

Thus the non-hermit icity of H is harmless. But in the exact expression (|5|), the hyperfine 
operator will remain in the orthogonality relations. 

One may also cast ([I]) into a strictly Hamiltonian form by introducing a secondary 8- 
component spinor ip s : 

(E - 2ap)V> = ips , Ei> s = (ml + 2K I )i). (45) 

This method is known from the relativistic treatment of spinless particles (for example from 
the Klein-Gordon equation). 

For V = —a/r, orthogonality relations are most elegantly derived in a dimensionless 
scaled variable, 

r = Er, d/dr = E~ 1 d/dr, p = p/E. (46) 
Dividing equation ([!]) by E 2 and setting E 2 = s for convenience, one obtains: 

V>(r) = 0. (47) 
Using the corresponding equation for x , one arrives at: 

is" 1 - sj 1 ) J xl m 2 + ^d 3 r = , s t = E 2 , (48) 

x\ m\ ^jd 3 r = 8ij . (49) 



1 i 

pa + -m 2 + /s - - + V(f) (1 - ip(a x cr 2 ) 



Remembering m? + = ml + m\ + 2m\mi{5, this is a simple generalization of the static limit 
m\j m 2 = 0. For positronium, the small components do not contribute to (|49f) (a previously 
proposed substitution r = Ep || gives more complicated orthogonality relations). 

Equation (|47| ) is explicitly CPT- invariant: Every bound state Sj has two different eigen- 
values Ei, namely Ei = y/si = m^i and Ei = —^fsl = — m^i, where m^i denotes the atomic 
mass in the state i (an excited atom is heavier than its ground state). The later value be- 
longs to the antiatom of mass m^i, i.e. m^i = rriAi- This follows from the CPT-invariance of 
QED, which ensures that the two-particle scattering amplitude at negative E describes the 
scattering amplitude of the two antiparticles. The range of the dimensionless radial variable 
f is < f < oo both for atoms and for antiatoms. In the old variable r, antiatoms have 
negative distances. This throws new light also on the static limit mi/m2 = 0. Here one 
defines E e = E — rri2 as the electron energy. For E < 0, one may use E e = E + m 2 . From the 
static Dirac equation in the variable f, one obtains a spectrum which is symmetric around 
E — ||, for V = —a/r. If one wants to keep this symmetry as a result of CPT also in the 
case of a finite nuclear charge distribution, one should parameterize V(r) rather than V(r). 

Of course, the mere CPT-invariance of a spectrum does not guarantee its correctness. 
Division of the Dirac-Breit equation by E and reformulation in terms of ?i = Ti/E, f 2 = 
r 2 /E, f = r/E also produces a CPT-invariant spectrum. But as the interaction in this 
case does not reproduce the QED Born approximation at all energies, one may hope that 
an 8-component formalism is again more successful. The corresponding equation has been 
presented in (|H]) and the Born approximation has been given in a suitable form in (^), but 
some details are still missing. However, it is clear that the 8-component formalism will be 
quadratic in the external potential V e , but linear in V. 
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APPENDIX A: COMBINATIONS OF THE HYPERFINE OPERATOR WITH ap 



Writing o^p as |(<r + Acr)p, one observes from ( pID that Act and cr x transform triplets 
Xt into the singlet Xs and vice versa. As a result, the combination required in ([5]) may be 
written as: 



^pA<r -iVpa x /E 



\ + {Pt- Ps)V/E 



pAcr. 



(Al) 



For total angular momentum /, the triplet states with I = f are excluded from pAcr by 
parity conservation. Thus one has in the notation of ||: 



pAcr = 2i 



f 


























\d + 


F/r 






\ 



p<T' 



( o o o -a_ ^ 

F/r 

V d + F/r / 



(A2) 



with p = rp/r , d± = d r ± 1/r and F — J f(f + 1). 
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APPENDIX B: CONNECTION WITH THE FORM DERIVED FROM BREIT 

OPERATORS 



We substitute in @ if} — e x ipB and multiply the equation by e~ x from the left, where 
the operator x is of the order of V/E and commutes with 75 and (3. To order a 4 , one may 
then approximate: 

e~ x pa e x « (1 - x) pa (1 + x) ^ pa + [pa, x] , e~ x Kj/E e x « Kj/E . (Bl) 
Choosing now 

x = -tr x tr 2 V/2E (B2) 
and extracting a common factor 75, one has: 

e~ x (pa, - iVpcr x /E) e x = p<n - ^(z<r x - o- 2 )[V, p] . (B3) 

The second piece is the hyperfine operator derived from Breit operators 0] and is known 
to reproduce the hyperfine structure of leptonium to order a 4 , including positronium ||. It 
can be rewritten in compact form: 

- -L(«r* +«r a )[V; V] = -L^V.V]^. (B4) 



APPENDIX C: SPIN SUMMATION 

The propagator of a lepton-antilepton pair will be needed in the perturbative interaction 
with radiation. Here we merely perform the spin summation for the trivial case of a free 
pair. We remind the reader that in the 16-component formalism, one defines for particles 
i = 1,2: 

7° = A , 7» = 7 4 V = ih^i , i>i = vai , (ci) 

which leads to the following form of the spin summation: 

J2 U l U 2UlU 2 = (fa + m 1 )(fa + m 2 ) . (C2) 

spins 

A similar notation may also be used in the 8-component version, but with the understanding 
7i = 7° = P 7 75i = 752 = 75 • Consequently, fa fa ^ fafa, but: 

faPfa = faPfa, fafa$ = Pfafa- (C3) 

The free leptonium equation (pT0|) and the corresponding equation for iy + become: 

{m 2 f3fa — mi^)w = , w + {m 2 f3fa — m^) = . (C4) 
From (|C3|) and (|C4|) , one easily verifies the following spin summation: 

s = ^2 uw ^ = m 2$i(3 + vn\fa . (C5) 



It is remarkable that this 8 x 8- matrix is linear in fa and fa, while the 16 x 16-matrix (|C2 ) 
also contains fafa. 
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